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Abstract: Random one plus two-body hamiltonians invariant with respect to 0(N\) © 
0(W 2 ) symmetry in the group-subgroup chains U(N) D U(M.) © U(N 2 ) 3 0{N\) © 0(W 2 ) 
and U(N) D 0(N) D O(Afi) © 0(./V 2 ) chains of a variety of interacting boson models are 
used to investigate the probability of occurrence of a given (wioj 2 ) irreducible represen- 
tation (irrep) to be the ground state in even-even nuclei; [oj\\ and [w 2 ] are symmetric 
irreps of 0(Ni) and 0(A/" 2 ) respectively. Employing a 500 member random matrix en- 
semble for N boson systems (with N = 10 — 25), it is found that for A/i,A r 2 > 3 the 
(cjiw 2 ) = (00) irrep occurs with ~ 50% and (wi^ 2 ) = (NO) and (07V) irreps each with 
~ 25% probability. Similarly, for A/i > 3,7V 2 = 1, for even N the w 1 — occurs with 
~ 75% and u\ = N with ~ 25% probability and for odd N, w x = occurs with ~ 50% and 
w\ = \,N each with ~ 25% probability. An extended Hartree-Bose mean-field analysis 
\ is used to explain all these results. 

Key words: Random interactions, nuclei, interacting boson models, group symmetries, 0(Afi)(B 
\ 0(A/" 2 ), ground state irreps, probabilities, mean-Held analysis. 

O ; 

^ ■ 1. INTRODUCTION 

Two-body random matrix ensembles (TBRE) defined over Hilbert spaces of various 
nuclear models led to the discovery that many of the regular features observed in low- 
lying levels and near the yrast line in nuclei can arise due to random interactions 
(with rotational symmetry) and this is opposed to the conventional ideas of using reg- 
ular (or coherent) interactions like pairing etc. in the nuclear hamiltonian. For the 
first time this result is found by Bertsch et al [1] using the shell model who showed 
that with random interactions ground states in even-even nuclei will be + with very 
high probability and they also generate odd-even staggering in binding energies, the 
seniority pairing gap etc. Similarly, Bijker and Frank [2] using the interacting boson 
model showed that random interactions generate vibrational and rotational structures 
with high probability. These unexpected results gave rise to a new field of research 
activity with random interactions in nuclei (they go beyond the TBRE applications for 
smoothed (with respect to energy) state densities, strength sums, transition matrix el- 
ements, information entropy in wavefunctions etc. in nuclei and other finite quantum 
systems; see [3] and references therein). In particular: (i) Zelevinsky and collabora- 
tors proposed the idea of geometric chaos for describing regular features generated 
by TBRE's; (ii) Arima's group introduced a variety of prescriptions, for predicting the 
probabilities, for simple systems such as single j-shell for fermions, single £-shell for 
bosons and some of their extensions; (iii) Bijker and Frank used a mean-field analy- 
sis with projective coherent states for interacting boson systems. For details of these 
studies we refer the readers to two recent reviews on this subject [4, 5]. 

1 Invited talk in the International symposium on 'Nuclear Physics' held at Dalian (China) 
during November 29 to December 1, 2003. Paper for the proceedings and the proceedings to 
appear as a special issue of High Energy Phys. & Nucl. Phys. (China). 
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Avery important aspect of TBRE's is that they admit group symmetries [6]. With m 
particles (fermions or bosons) in Af single particle states there is a U (77) spectrum gen- 
erating algebra (SGA). In all the shell model/Interacting boson model studies reviewed 
in [4, 5] one plus two-body hamiltonians that are rotational scalars are considered, i.e. 
all the 0(3) scalars in U(Af) D 0(3) with one and two particle matrix elements of the 
one and two-body parts respectively chosen to be random variables (in some studies 
the one body part is dropped). Immediately one sees that TBRE's can be extended 
to scalars of various subgroups of (7(77) > i.e. scalars of G in (7(77) D G D . . . D 0(3). 
The purpose of the present paper is to consider such an extension to 0(N\) © 0(7V 2 ), 
M + Mi = M which appears in a very large class of interacting boson models (IBM's) 
used in nuclear structure and address the question of with what probability a given 
(u>icj 2 ) irreducible representation (irrep) will be the ground state in even-even nuclei; 
note that [wi] and [co 2 ] are symmetric irreps of 0(M) and 0(77 2 ) respectively. In Section 
2 given are the random one plus two-body hamiltonians with 0(N\) © 0(J\f 2 ) symmetry 
in interacting boson models. Section 3 gives the results of numerical TBRE calcula- 
tions and their understanding using an extended Hartree-Bose mean-field analysis. 
Finally Section 4 gives conclusions and future outlook. 

2. RANDOM INTERACTIONS WITH 0(M) © 0(7V 2 ) SYMMETRY IN IBM's 

Large class of interacting boson models (IBM's) of nuclei admit U(Af) D U{N\) © 
U{N 2 ) D 0(M) © 0{N 2 ) and [/(TV) D 0{M) D 0(7Vi) © 0{M 2 ) group-subgroup chains; 
J\f = Wi +77 2 . Examples (all for even-even nuclei) are: (i) spIBM or nuclear vibron model 
[7] with C/(4) SGA and (77i,77 2 ) = (3,1); (ii) sdlBM for quadrupole collective states [8] 
with U(6) SGA and (7Vi,7V 2 ) = (5, 1); (iii) spdlBM for GDR states [9] with {7(9) SGA and 
(7Vi,7V 2 ) = (8,1), (6,3), (5,4); (iv) sdglBM for quadrupole plus hexadecupole states [10] 
with 17(15) SGA and = (14, 1), (9, 6), (10, 5); (v) sdpflBM for octupole states [11] 

with 17(16) SGA and (M,A/" 2 ) = (15,1), (10,6) etc; (vi) sdgpflBM [12] with C/(25) SGA 
and (7Vi,7V 2 ) = (24,1), (15,10) etc; (vii) spp'IBM or the U(7) model for 3-body clusters 
in nuclei [13] with U{7) SGA and (7Vi,7V 2 ) = (6, 1), (4,3); (viii) IBM-3 or the isospin (T) 
invariant seflBM (here the bosons carry T = 1 degree of freedom) [14] with J7(18) SGA 
and (77i,7V 2 ) = (15,3); (ix) IBM-4 or the spin-isospin [S,T] invariant seflBM (here the 
bosons carry (ST) = (10) © (01) degree of freedom) [15] with (7(36) SGA giving examples 
with {Afi,Af 2 ) = (30,6), (3,3), (18, 18), (15, 15); (x) IBM-2 or proton-neutron IBM [16] with 
U(12) SGA and (Wi, Af 2 ) = (10, 2). In this paper, for simplicity, we consider group chains 
with Wi > 3,7V 2 = 1 and Wi > 3,W 2 > 3, i.e. 7V 2 = 2 situations (as in (x) above) are not 
considered. 

Group chains, for symmetric U(Af) irreps {N} the irrep labels for other group algebras 
in the chains and their reductions for the Afi > 3 and M 2 = 1 situation (hereafter called 
I) are, 



In (1,2), label(s) a for the irreps of K need not be specified as the algebra K do not play 
any role in the present work. Note that (7(77) D f7(M = TV - 1) © U{N 2 = 1) and the 




and 




(2) 
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U(J\f 2 = 1) and its irreps {n 2 }, n 2 = N - m are not shown in (1). The general one plus 
two-body 0(M) (0(7V 2 = 1) will not exist) scalar hamiltonian built out of the Casimir 
operators of the group algebras in the chains (1,2) is given by 

H 1 = ±[a 1 C 1 (U(M 1 ))+a 2 C 1 (U(Af 2 = l))} 

+WW=Tj [a 3 C 2 (?7(M)) + a 4 C 2 (U(Af 2 = 1)) + a 5 Ci([/(AA 1 ))C 1 (C/(AA 2 = 1)) (3) 
+ a 6 C 2 (0(Af)) + u 7 C 2 {0{Mi))] 

Similarly for the J\f\ > 3 and Af 2 > 3 situation (hereafter called II), the group chains, for 
symmetric U(Af) irreps {TV} the irrep labels for other group algebras in the chains and 
their reductions are, 

U{M) D U(Ni) e U{N 2 ) D 0(M) © 0{N 2 ) D K \ 
{N} {m} {n 2 } [ Wl ] [uj 2 ] a I 

m = 0, 1, 2, • • • , N; n 2 = n-nx 

u>\ = ni,ni — 2, . . . ,0 or 1, uj 2 = n 2 , n 2 — 2, . . . , or 1 

and 

U{N) D 0(Af) D 0(M) © 0(Af 2 ) D K \ 
{iV} [w] [wi] [wa] a / 

w = iV,iV-2,...,0orl, wi +w 2 = w,w - 2, . . . ,0 or 1 (5) 

The general one plus two-body 0(N\) 0(Af 2 ) scalar hamiltonian built out of the 
Casimir operators of the group algebras in the chains (4,5) is given by 

H" = i[a 1 C 1 ([/(AA 1 )) + a 2 C 1 (C/(AA 2 ))] 

+jv07=ij ["3C 2 (f/(M)) + a±C 2 {U{M 2 )) + o B Ci(i;(M))Ci(^(M)) (6) 
+ a 6 C 2 (0(AA)) + a 7 C 2 (0(M)) + a 8 C 2 (0(AA 2 ))] 

In (3,6), Ci and C 2 are linear and quadratic Casimir invariants and their matrix el- 
ements for example are (Ci(C/(M))) {ni} = m, (C 2 {U - ni(m +M - 1) and 
(C 2 (0(7Vi))) kl1 = wi(wi +M - 2). Given N bosons, in the \Nrnui) basis for I and 
\N(nin 2 )(u)ioj 2 )) basis for II, the many boson H matrix for H 1 and H n respectively will 
be always tridiagonal; the a 6 terms in Eqs. (3,6) generate off-diagonal matrix elements. 
In particular I is a generalization of the splBM analyzed before [17, 18] while II is com- 
pletely new. For each allowed ui in I and (wiw 2 ) in II the H matrices are constructed 
using the transformation brackets, given in [19], between the chains (1) and (2) for I 
and (4) and (5) for II. The matrices are diagonalized, for each member of a 500 member 
TBRE, for boson numbers N = 10 — 25. Thus, in the calculations the parameters in 
Eqs. (3,6) are chosen to be independent Gaussian variables with zero mean and unit 
variance and 500 samples of the same are considered. Now we will discus the results. 

3. RESULTS FORTVi > 3,W 2 = 1 AND M > 2>,N 2 > 3 SYSTEMS AND THEIR 

MEAN-FIELD ANALYSIS 

3. 1 Results of Numerical Calculations 

Figures la and lb give the probabilities, in the situation A/i > 3 and M 2 = 1, i.e. for I 
as M\ is varied. In Fig. la shown are the results for ui = 0, N for even boson number 
N [N = 10) and in Fig. lb for wi =0,1,N for odd N (N = 15) to be ground states. Fig. 
lc shows the same results but as a function of the boson number N for A/i = 14. It 
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should be noted that the probabilities are negligiblly small for the [loi] irreps not shown 
in the figures. In general for even N, uj\ = is ground state irrep with ~ 65 - 74% and 
lui = N with ~ 25 - 32% probability. Similarly for odd N, u)\ = is ground state irrep 
with - 46 - 50%, wi = 1 with - 18 - 24% and wi = N with - 24 - 30% probability. They 
reproduce the spIBM results known before [17, 18, 5] and provide a test of the present 
calculations. They are also close to the seflBM results known before (although in these 
studies K = 0(3) is chosen and H 1 to be a 0(3) scalar, note that wi = 1 gives L = 1 in 
splBM and L = 2 in sdlBM) [18, 20, 5]. 
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Fig. 1. Probabilities for various group Irreps to be ground states, (a) and (b) give the variation 
as a function of M for fixed boson number and (c) and (d) as a function of the boson number N. 
(a), (b) and (c) are for I and (d) is for II. See text for details. 

In the situation Wi > 3, W 2 > 3 (i.e. II), Figure Id (for various TV values with (7Vi,7V 2 ) = 
(9,6)) and Table 1 (for various (A/i,W 2 ) with N = 10,20) give the results for the irreps 
(wiw 2 ) = (00), (M)), (ON). Here only even N is considered. It is seen that in general 
(wiw 2 ) = (00) is ground state with — 50 - 55% and (wiw 2 ) = (ON) and (NO) each with 
<~ 20 — 24% probability. The probabilities for other (wiw 2 ) irreps to be ground states is 
negligiblly small. Before giving a mean-field analysis of these results some remarks 
are in order. In sdglBM with (M,^) = (9,6), for 7-soft nuclei the (wiw 2 ) = (ON) is 
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expected to be the ground state but, as seen from Fig. Id, random interactions give 
this irrep only with ~ 20% probability. Therefore random interactions are not good for 
sdglBM for (M,M 2 ) = (9,6). However for (Mi,M 2 ) = (14, 1) (this system is used recently 
in phase transition studies [21]), the U(1A) irrep [0] occurs, for even N, as ground state 
with ~ 70% and thus random interactions may be useful here. Another example is, in 
IBM-4 with (Mi,M 2 ) = (3,3) the ( Wl w 2 ) = (ST) = (00) irrep is ground state with - 50% 
probability. However in real nuclei this irrep is expected to be the ground state and 
therefore in IBM-4 one can use random interactions but a regular part enhancing the 
probability for (ST) = (00) should be added. 

3.2 Mean-field analysis 

Bijker and Frank [18, 20] carried out a mean-field analysis of sp and sd IBM's to give 
a quantitative understanding of the probabilities, with random interactions in these 
models, for a given L to be the ground state. We will follow this approach with suitable 
extensions to describe the results found in Fig. 1 and Table. 1. We begin with I, i.e. 
Mi > 3, M 2 = 1. The M 2 = 1 gives s bosons (with angular momentum I = 0). Similarly 
M > 3 gives bosons carrying M\ degrees of freedom and they can be thought of as 
bosons with I = £\,l 2 , ■ ■ -i r > Z)[=i( 2 ^ + 1) = M- Just as in [20], a one parameter H is 
considered (with a\ and a e terms in (3)), 

H = ± cosxfi! + x) sinx S+S- , S + = s V - £ *>{ ' b{ , S. = (S+)t (7) 

In (7) — 7r/2 < x < 37r/2 so that all attractive and repulsive interactions are included. 
The ground state shapes and hence the ground state uji are determined by minimizing 
the energy functional for the axially symmetric coherent state (CS), 

\N a) — -L= ^cos a + sin a x^j |0) , xj = -^^6j. (8) 

In Eq. (8), -n/2 < a < ir/2 and it gives for example correctly the CS used for sdglBM in 
the past [10]. The energy functional E(a) = (N a \ H \ N a) = cos\ sin 2 a+j sinx cos 2 2a. 
The minima of E(a) and the corresponding shape parameters a of the CS divide into 
three classes: (i) a = for -n/2 < x < t/4, i.e. in a 37r/4 range; (ii) a = ir/2 for 37r/4 < 
X < 3tt/2, i.e. in a 37r/4 range; (iii) a such that cos 2a = cotx for 7r/4 < x < 37r/4, i.e. in 
a 7r/2 range. It is seen from (8) that a = gives s boson condensate and hence here 
loi = and it occurs with (37r/4(27r)) x 100% = 37.5% probability. Similarly a = ir/2 gives 
x boson condensate and here, apart from a constant factor, E = — sin x ^i (t<->i +Mi — 2) 
where uj\ = m,n\ — 2, ... ,0 or 1. Then clearly, sinx positive (this happens 7r/4 times) 
gives uj\ = N to be lowest with 12.5% probability and sinx negative (this happens n/2 
times) gives uj\ = to be lowest for even and w\ = 1 for odd N with 25% probability. 
For cos 2a = cot x the condensate is deformed with both s and x bosons. It is plausible 
to argue that the condensate here gives a band with L = kwi, u>i = 0, 1, 2, . . . , TV and 
k = 1 for spIBM, k = 2 for srflBM, k = 4 for sdglBM etc. The moment of inertia (T) of 
these bands should follow from 0(M\) cranking (a method for this may be possible via 
the results in [22]). As yet there is no theory for this and therefore we assume that the 
0(3) cranking formula given in [18, 20] is valid here to within a constant. Then (with 
T = (sin x — cos x) / sin \ cos x) it is easily seen that for cos 2a = cot x, u\ = is lowest with 
12.5% probability and uji = N is lowest with 12.5% probability. Combining all these 
results will give for I for the ground state probabilities: (i) ui = with 75% and uj\ = N 
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with 25% for even N; (ii) uj x = with 50%, wi = 1 with 25% and ui\ = N with 25% for odd 
N. They give a good description of the results in Figs, la, lb and lc. 



Table 1. Probabilities (in percentage) for (wiw 2 ) to be ground state irrep 



Model 


M 


A/2 


N 


(WHU2) = (00) 


(u>iw 2 ) = {NO) 


(W1W2) = {ON) 


U{7) 


4 


3 


10 


55.4 


21.4 


20.5 








20 


54 


21.2 


20.6 


spdYBM 


6 


3 


10 


55 


22.4 


19.8 








20 


53.5 


21.9 


20.3 


sdgYBM 


10 


5 


10 


55.3 


22.9 


19.2 








20 


53.6 


22.8 


20.4 


spdflBM 


10 


6 


10 


49.3 


24.6 


21.9 








20 


49.3 


23.8 


22 


sdgpflBM 


15 


10 


10 


53.8 


22.9 


20.3 








20 


54.4 


22.4 


20.5 


IBM-3 


15 


3 


10 


49 


27.1 


19.8 








20 


49 


25.6 


20.7 


IBM-4 


3 


3 


10 


49.2 


22.8 


22.8 








20 


48.8 


23.2 


23.2 




30 


6 


10 


50.1 


28.6 


18.6 




18 


18 


10 


50 


23.5 


23.5 




15 


15 


10 


49.6 


23.6 


23.6 



Now we will consider the mean-field analysis for II, i.e. for A/i, A/2 > 3 and the discus- 
sion will be restricted to even N. Just as the operator in Eq. (8), let us introduce 
yt and *t operators, y\ = ^ £Li b\ ifi , £? =1 (24 + 1) = M and 4 = _I_ £* =i ^ 

Ej=i( 2 ^- + !) = A/2. Then the hamiltonian, CS and E{a) are, 
= ± cosx" 2 + jy^zT) sinx S+S- , 

5+ = 5+(l) - S+(2) = ^ 6j ( • 6j ( 6], , S_ = (5+)t 

fci 3 3 (9) 

\ N a ) = ^m( cosa yo + sinaz o) 1°) 

-E(a) = cosx sin 2 a + i sinx cos 2 2a 

The i? in (9) is same as H 11 defined in Eq. (6) but only with the a 2 and a 6 terms. 
The equilibrium shapes correspond to a = 0, a = ir/2 and cos 2a = cotx with the range 
of x' s j us t as before. The a = gives y-boson condensate with energy E(a = 0) cx 
— sinxwi(wi +A/i - 2). Then the ground state irreps are (wiw 2 ) = (00) with 25% and 
(W1W2) = (A/0) with 12.5% probability. Similarly a = 7r/2 gives z-boson condensate with 
energy E(a = tt/2) cx - sinx ^2 (^2+ A/2 -2) and then the ground state irreps are (wiw 2 ) = 
(00) with 25% and (wiw 2 ) = (OA/) with 12.5% probability. In the situation cos 2a = cotx, 
cranking has to be done with respect to both 0(N\) and 0(A/2). Evaluating moment of 
inertias as before gives E to be, to within a constant, 

o>i(ui +M - 2) o) 2 (^2 +A/" 2 - 2) (sin x± cosx) n nl 

^ = 1 -. , A ± = T : U<->J 

A + cosx sinx 

With 7r/4 < x < 37r/4 here, it is seen that A + is +ve and A_ is — ve for 7r/4 < x < t/2 
and A + is -ve and A- is +ve for tt/2 < \ < 37r/4. Therefore, here (A/0) and (OA/) irreps 
will be ground states each with 12.5% probability. Combining all the results give for 
II, (wiw 2 ) = (00), (ATO) and (OA/) irreps to be ground states with 50%, 25% and 25% 
probability. These numbers clearly describe the results in Table 1 and Fig. Id. 
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In summary, the mean-field approach of [18, 20] with proper extensions gives a good 
understanding of the results in Fig. 1 and Table 1 although all the results are obtained 
using a constant probability for \ in Eqs. (7,9). Extension of the analysis to odd N for 
II (here the lowest (wiw 2 ) are (10) and (01)) and also the calculations for Af 2 = 2 will be 
given elsewhere. 

4. CONCLUSIONS AND FUTURE OUTLOOK 

In this paper for the first time TBRE's preserving irreps of group symmetries (other 
than 0(3)), for boson systems, are introduced and showed that the + dominance 
observed in ground states extends to group irreps. An extended mean-field analysis 
is shown to give good description of the numerical results obtained for a variety of 
interacting boson models. The mean-field analysis in Section 3 is restricted to the 
simple mixing hamiltonians given by Eqs. (7,9) and in a future paper this will be 
extended to the full hamiltonians given by Eqs. (3,6). Similarly, in a future publication 
we will consider 0(N), N > 3 cranking so that the application of 0(3) cranking to I and 
II in Section 3 can be validated. At present the justification for using (9(3) cranking 
comes from the good agreement between the mean-field results and those in Fig. 1 and 
Table 1. It should be added that Kusnezov's analysis for splBM [17], based on random 
polynomials, can be applied to H 1 and H 11 as the matrices here are tridiagonal and 
this will be done elsewhere. Finally, it will be interesting to extend the present work to 
other general classes of group-subgroup chains in IBM's (see [14, 15] for examples) and 
also to group chains for fermion systems (as they appear for example in the shell model 
[23]). It is plausible that the results of these extensions will give deeper understanding 
of geometric chaos and regularities generated by random interactions. 

A. Arima's talk in 'Symmetries in Science XIW held at Bregenz in July 2003 has pro- 
vided the basis for initiating the work presented in this paper. The author thanks N. 
Yoshinagafor discussions in the Bregenz meeting and Y.M. Zhao for very useful corre- 
spondence. This article is dedicated to Prof. J.P. Draayer on his 60 th birthday. 
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